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The statistical properties of turbulence are considered to be universal at sufficiently small length 
scales, i. e., independent of boundary conditions and large-scale forces acting on the fluid. Analyzing 
data from numerical simulations of supersonic turbulent flow driven by external forcing, we demon- 
strate that this is not generally true for the two-point velocity statistics of compressible turbulence. 
However, a reformulation of the refined similarity hypothesis in terms of the mass- weighted velocity 
p^^^v yields scaling laws that are almost insensitive to the forcing. The results imply that the most 
intermittent dissipative structures are shocks closely following the scaling of Burgers turbulence. 
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The notion of universality is central to theoretical and 
observational accounts of turbulence. For incompress- 
ible turbulence, the mathematical analysis carried out 
by Kolmogorov led to the famous 2/3 law for the second 
order structure function of turbulent velocity fluctuations 
P, m . Remarkably, this scaling law was experimentally 
confirmed even if the premises of Kolmogorov's theory- 
statistical equilibrium, homogeneity and isotropy-were 
not satisfied. It is commonly accepted that this is due 
to an inertial subrange of scales, where the dynamics of 
turbulence conditions the flow such that these premises 
are met asymptotically toward smaller length scales inde- 
pendent of the large-scale properties of the flow. Whereas 
terrestrial applications are mostly concerned with the in- 
compressible regime, the idea that universal scaling also 
exists in the inertial subrange of the highly compress- 
ible, supersonic turbulence has become popular in as- 
trophysics. Particularly, observational properties of star- 
forming clouds are explained by supersonic turbulent mo- 
tion which is seen as an agent that controls the forma- 
tion of stars besides gravity [3]. We tested the hypothe- 
sis of universality on data from numerical simulations of 
supersonic turbulence with the same root mean square 
(RMS) Mach numbers but different large-scale forcing. 
The result is that universality in the sense of the Kol- 
mogorov theory is clearly violated even if intermittency 
corrections are applied, but scaling laws that are nearly 
indepenent of the forcing apply to the mass- weighted ve- 
locity V = p^l'^v introduced by Kritsuk et al. [4]. More- 
over, the scaling exponents turn out to be consistent with 
log-Poisson models, if the most intermittent dissipative 
structures are shocks fulfilling the scaling of Burgers tur- 
bulence. 
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The fundamental relation for the scaling properties of 
incompressible turbulence is the refined similarity hy- 
pothesis P, Q , 

S^{l)=C^P'l\e^l^), (1) 

where the are constant dimensionless coefficients, 
Srp(P) := {^v\) are the statistical moments of the ve- 
locity fluctuation = \v{x^t) — v{x -\- i^t)\ of order 
p = 1,2,3,... and is the rate of energy dissipation per 
unit mass averaged over a region of size i. The brackets 
( ) denote the ensemble average. The factor (ef ^) ex 
is attributed to the intermittency of turbulence [1]. The 
scaling exponents of Sp{i) in the inertial subrange are 
thus given by Cp = p/3 + rp/3. 

Modeling turbulent energy dissipation by a ran- 
dom cascade obeying log-Poisson statistics, Dubrulle 
showed that the relative scaling exponents Zp have the 
general form 

Z,,= | = (l-A)| + :^(l-;3^'=)^ (2) 

The intermittency parameter P is interpreted as a ran- 
dom cascade factor relating dissipative structures of dif- 
ferent intensity, and C = A/(l — /?) is the co-dimension 
of the most intense dissipative structures [?(. The co- 
dimension is related to the fractal dimension by C = 
1^ — 3. As argued by She and Leveque the scaling of 
these structures, £~^^ is given by the inverse of the ki- 
netic energy available for dissipation at the length scale 
£. For incompressible turbulence, the most intense dissi- 
pative structures are assumed to be vortex filaments, for 
which C = 2 and A = 2/3. 

Boldyrev joj] proposed C = 1 for supersonic turbulence 
(Kolmogorov-Burgers model), because he considered the 
most intense dissipative structures to be shocks, while 
keeping A = 2/3 as in the She-Leveque model for in- 
compressible turbulence. However, the kinetic energy at 
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FIG. 1: Interpolated relations between the time- averaged 
structure functions Sp and Ss for solenoidal (solid lines) and 
compressive (dashed lines) forcing. The order p ranges form 
1 (bottom) to 5 (top). The corresponding power- law fit func- 
tions are plotted as thin lines. 



the length scale i is proportional to i for Burgers tur- 
bulence {Sv£ (X Following the arguments by She 
and Leveque we propose that the most intense dissipa- 
tive structures should obey the scaling law rather 
than i. e., A = 1. The scaling exponents obtained 
from equation ([2]) for A = 1 are markedly different from 
the prediction of the Kolmogorov-Burgers model. 

In the following, we will determine the relative scal- 
ing exponents := Cp/Cs from the relations S^(P) = 
Ss{£)£^p fld\ for simulations of supersonic isothermal tur- 
bulence with periodic large-scale stochastic forcing [11^ 
In these simulations, the compressible Euler equations 
were solved with the piecewise parabolic method in the 
FLASH3 implementation jTH, E^. In one case, purely 
solenoidal (divergence-free) stochastic forcing was ap- 
plied, in the other case the forcing was purely compres- 
sive (rotation- free) [13, [lH . We use the term compres- 
sive synonymous to dilatational (rotation- free). In each 
case, the system was evolved over 10 auto-correlation 
time scales (subsequently denoted by T) of the force field 
at grid resolutions TV = 256^, 512^ and 1024^. Turbu- 
lence was found to be fully developed with a steady-state 
RMS Mach number 5.5 after about two autocorrela- 
tion time scales. Following previous numerical studies 
of supersonic turbulence [16, 17], we computed transver- 
sal structure functions Sp{£), i. e., the p-th moments of 
the velocity fluctuation projected perpendicular to the 
spatial separation £. The structure functions were com- 
puted from a statistically converged sample in the inter- 
val 2 < t/T < 10 using a Monte Carlo algorithm [IL] . 



TABLE I: Relative scaling exponents Zp from fits of time- 
averaged structure functions Sp vs. . 





N 


Zi 


Z2 


Zs 


Z4 


Zs 




256^ 


0.472 


0.786 


1. 


1.160 


1.289 


sol. 


512^ 


0.474 


0.792 


1. 


1.149 


1.265 




1024^ 


0.466 


0.788 


1. 


1.150 


1.266 




256^ 


0.603 


0.879 


1. 


1.072 


1.126 


comp. 


512^ 


0.627 


0.896 


1. 


1.056 


1.097 




1024^ 


0.628 


0.897 


1. 


1.055 


1.095 



The structure functions Sp{i) averaged over the time 
interval 2 < t/T < 10 for TV = 1024^ are plotted as 
functions of the time-averaged third-order structure func- 
tion S^{i) in Fig. [H The exponents Zp are given by 
the slope of logS^ vs. logS'g^. For the determination 
of linear fit functions, fitp ex ZplogS^^ we imposed the 
criterion Vp < 5 : err^ := | exp (fit^) - S^\/S^ < 0.01 
in the fit range. This error criterion is fulfilled in the 
intervals 12.0 < 5*3^ < 120 for solenoidal forcing and 
25.0 < 6'3^ < 150 for compressive forcing. One should 
note that the relations between Sp{£) and S^{£) agree 
quite closely with the fit functions (err^ < 0.05) for large 
length scales corresponding to the maximal values of . 

The relative scaling exponents Zp inferred from the 
time-averaged transversal structure functions are listed 
in Table [D In the case N = 1024^, the standard errors 
of the parameters Zp are of the order 10~^. We estimate 
systematic errors to be of the order 10~^. In Fig. [21 we 
compare these values to the instantaneous scaling expo- 
nents as functions of time for N = 1024^. Most impor- 
tantly, we see that the scaling exponents resulting from 
solenoidal forcing differ markedly from the case of com- 
pressive forcing. In each case, the variation of Zi and Z2 
over the time scale T clearly shows temporal correlation 
and there appears to be anti-correlation with the scal- 
ing exponents of order greater than three. This suggests 
that the instantaneous scaling exponents show an im- 
print of the stochastic variation of the large scale forcing 
rather than purely statistical scatter. Even for incom- 
pressible turbulence, an influence of the large scales on 
much smaller scales is reported [18]. 

The Kolmogorov-Burgers model implied by equa- 
tion ([2]) for A = 2/3 and /3 = 1 - A/C = 1/3 is 
plotted as dotted line together with our data in Fig. [3l 
Clearly, there are large deviations both for solenoidal 
and for compressive forcing. One reason is that the 
Kolmogorov-Burgers model only applies in the hyper- 
sonic limit, whereas the most intense dissipative struc- 
tures are constituted by varying fractions of vortex fil- 
aments and shocks depending on the RMS Mach num- 
ber. Then one would expect 1 < C < 2 [16]. More- 
over, fitting the log-Poisson model with A = 2/3, we 
find Csoi ^0.76 and Ccomp ^ 0.67, which is about the 
minimal co-dimension C = A = 2/3, for which the Zp are 
real. In the case of compressive forcing, no closely match- 
ing fit function exists. On the other hand, fitting the 
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FIG. 2: Time evolution of the relative scaling exponents Zp. The time averages are indicated by the horizontal lines. 





FIG. 3: Comparison of the relative scaling exponents Zp 
for solenoidal and compressive forcing (see Table |l| with 
the Kolmogorov- Burgers model (dotted line) and log-Poisson 
model fits with A = 2/3 (dot-dashed lines), C — 1 (thick 
dashed lines) and A = 1 (solid lines). The vertical bars indi- 
cate the standard deviations of the instantaneous values. 



one-parameter family of models with A = 1, we obtain 
co-dimensions Cgoi ^1.5 and Ccomp ^1.1. These mod- 
els match the time-averaged relative scaling exponents 
very well. We cannot fully discriminate other families of 
models though. For instance, assuming C = 1 as in the 
Kolmogorov-Burgers model and varying A as fit param- 
eter, yields Agoi = 0.79 and Acomp ^ 0.94. Nevertheless, 
in the case of compressive forcing, the match with the 
numerically computed scaling exponents is much better 
in comparison to the A = 2/3 models. 

To carry over the log-Poisson models to compressible 
turbulence, the refined similarity hypothesis might be 
applied in the form ([T]), where the rate of energy dis- 
sipation at length scale i is understood to be the Favre- 
filtered quantity ei — {pe)i/ p^. The filter operation ( )i 



FIG. 4: Comparison of the relative scaling exponents Zp 
(mass- weighing v — p^^^v) for solenoidal and compressive 
for cine: fsee Table ini with various models analogous to Fig. El 



smoothes out fluctuations at length scales ^ ^, ^= {p)i 
and e{x^t) is the local rate of energy dissipation per unit 
mass. On the other hand, one might consider the rate 
of dissipation per unit volume, q := p^e^, as the vari- 
able from which the hierarchy should be constructed [l^ . 
This is also suggested by the alternative formulation of 
the log-Poisson model by She and Waymire • Following 
this proposition, the refined similarity hypothesis with 
mass- weighing. 



ss) = Cpt^'^'Hipt^tY'^)^ 



(3) 



where Sp{€) := {S{p^^^v)^)^ leads to the relative scalings 
Sp{i) = Ss{i)£^p for any (s > 0. As proposed by Kritsuk 
et al. [4], the exponents Zp associated with the two-point 
statistics of i; := p^/^v, are then given by an expression 
analogous to ([2|) with parameters /3 and A. 

The mass- weighted structure functions S^{£) vs. 
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FIG. 5: Relations between the time- averaged structure func- 
tions Sp and Ss with mass weighing analogous to Fig. [T] 



TABLE II: Relative scaling exponents Zp from fits of time- 
averaged mass- weighted structure functions Sp vs. with 



mass wei^ 


^hing V 














N 


Zi 


Z2 


Z3 


Z4 


Zs 




256^ 


0.546 


0.839 




1.094 


1.150 


sol. 


512^ 


0.550 


0.845 




1.082 


1.122 




1024^ 


0.539 


0.840 




1.080 


1.112 




256^ 


0.635 


0.893 




1.034 


1.026 


comp. 


512^ 


0.634 


0.887 




1.050 


1.068 




1024^ 


0.605 


0.869 




1.066 


1.100 



S^{i) computed from our numerical data are plotted in 
Fig. [5l In order to satisfy the chosen error criterion, 
we had to restrict the fit functions to narrower ranges 
15.0 < 5'3^ < 100 and 25.0 < ^'3^ < 75.0 for solenoidal 
and compressive forcing, respectively. The resulting val- 
ues of the relative scaling exponents Zp are summarized 
in Table nil Although the sensitivity on resolution is more 
pronounced than for the Zp and the scatter of the instan- 
taneous values of Zp is very large due to the pronounced 
intermittency of the mass density [17], trends are nev- 
ertheless discernable. Most noticeably, the differences 
between turbulence driven by solenoidal and compres- 
sive forcing are substantially reduced, especially, for the 
higher-order exponents. As shown in Fig.HJ again neither 
the Kolmogorov-Burgers model nor log-Poisson models 
with A = 2/3 match the scaling exponents Zp. The clos- 
est fits are obtained for A = 1, where Csoi ^ 1.18 and 
C'comp ~ 1.08, while the models with fixed co-dimension 
C = 1 yield Asoi = 0.90 and Acomp ^ 0.92. 

In summary, the scaling of turbulent supersonic veloc- 
ity fields is characterized by power laws that vary sub- 
stantially with the large-scale forcing. This finding dis- 
proves the conventional notion of universality in the su- 
personic regime and bears consequences on the theory of 
turbulence-regulated star formation. On the other hand, 
calculating two-point statistics of the mass-weighted ve- 
locity V = p^^^v, we found that the influence of the 
forcing was considerably reduced. Based on the corre- 
sponding formulation of the refined similarity hypoth- 
esis (j3j), the scaling exponents are very well described 
by log-Poisson models ([2]), for which the parameter A is 
close to unity. This is the expected value for the most 
intense dissipative structures in supersonic turbulence. 
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